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>  ^oSTRACT  ^  _  ,  ,  , 

i  Operational  and  tactical  calculations,  as  well  as  scientific 

}  and  technical  investigations,  often  involve  consideration  of  random 
phenomena  and  processes  which  can  he  characterized  by  statistical 
[  parameters  calculated  as  averages  of  observed  (selected)  values.  The 
■  statistical  probability  modelling  method  for  random  processes  suggested 
j  uses  the  Monte  Carlo  method  as  its  mathematical  base,  and  the  article 
proceeds  to  describe  its  application  to  two  practical  examples,  one  of 
*  which  is  the  calculation  of  the  process  involved  in  a  missile  closing 
a  target;  the  other  reproducing  an  actual  noise  process  in  accordance 
|  with  the  normal  distribution  law.  The  mathematical  method  described 
i  involves  the  use  of  an  electronic  computer  to  obtain  the  statistical 
j  evaluation  of  characteristics  with  required  accuracy  and  to  reproduce 
j  the  probability  models  of  the  subjects  under  investigation  using 
experimental  data  based  on  distribution  laws  and  correlated  ties  in  the 
elements  of  those  subjects,  without  waiting  for  the  development  of  an 
analytical  theory  and  the  disclosure  of  the  physical  content  of  those 
subjects.  Orig.  art.  has:  8  formulas  and  1  figure. 

/v 
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STATISTICAL-PROBABLE  MODELING  OP  RANDOM  PROCESSES 


A.  S.  Zhernenko 

During  operational-tactical  calculations,  as  well  as  during  the 
fulfillment  of  scientific-technical  investigations  it  is  often 
necessary  to  discuss  random  phenomena  and  occurrences. 

Parameters  of  the  movement  of  rockets  and  missiles,  characteristics 
of  a  naval  medium,  elements  of  motion  of  a  target  —  all  these  are 
values,  which  during  observations  (measurements)  may  acquire  various 
discrete  or  continuous  values. 

In  a  random  occurrence  instead  of  the  rigid  cause  of  communication 
with  a  fixed  complex  of  external  conditions  of  its  realization  is  only 
observed  a  statistical  bond,  which  is  revealed  at  a  greater  series  of 
observations  and  is  expressed  in  a  probable  form.  Correspondingly 
characteristic  feature  of  a  random  process  appears  in  the  fact,  that 
it  cannot  be  previously  predicted,  which  of  the  multitude  describing 
its  functions  will  be  observed  at  an  alternate  experiment. 

Random  values  and  processes  are  characterized  by  statistical 
parameters,  which  are  calculated  as  an  average  of  observed  (selective) 
values . 

A  random  value  £  is  considered  known,  when  all  the  values  are 
known,  values  which  it  can  accept  (x^,  x2>  ...»  xn...),  as  well  as 
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probabilities  of  realizing  each  of  these  values  (P  (v,).  p  P  (.v„)...). 

The  exhaustive  characteristic  of  a  random  value  £  serves  either  the 
function  of  probability  distribution  F(x),  or  the  density  of 
probabilities  W(x). 

An  arbitrary  function  (x(t)}  can  be  represented  as  a  sequence  of 
an  infinite  multitude  of  random  values  x(tk)  —  instantaneous  values 
of  its  realization  in  a  discrete  series  of  points  of  argument  (t). 
Hence  it  follows,  that  statistical  properties  of  a  random  process  can 
be  fully  determined  only  by  an  infinite  sequence  of  distribution 
functions.  To  obtain  the  mentioned  sequence  is  in  a  majority  of  cases 
practically  impossible.  Therefore  searches  are  conducted  for  methods 
of  investigating  random  processes,  allowing  to  confine  oneself  only 
to  investigations  of  distribution  functions  of  the  first  order.  The 
method  of  statistical-probable  modeling  of  random  processes,1 
pertaining  to  their  number,  is  based  in  its  mathematical  basis  on  the 
known  Monte-Carlo  method. 

A  statistical-probability  model  of  a  stationary  random  process 
is  represented  in  form  of  an  ensemble  of  discrete  realizations, 
obtained  on  the  basis  of  two  statistical  characteristics  of  the 
process:  correlation  function  (x)  and  monodimensional  distribution 

function  F^(x). 

The  correlation  function  as  a  time  averaged  statistical  parameter 
of  stationary  processes  is  equal  to  mathematical  exp'  ctation  of  a 
product  of  two  instantaneous  values  —  x(t)  and  x(t  +  x),  —  shifted 
in  time  to  interval  x.  Being  a  statistical  average  of  second  order, 
the  correlation  function  represents  a  simplified  reflection  of  time 


^hls  method  was  developed  at  the  Cybernetics  Inst,  of  the  Acad, 
of  Sc.  Georgian  SSR  under  the  leadership  of  Dr.  of  Phys.  Math.  Sc. 

V.  V.  Chavchanidze  and  explained  in  detail  in  his  reports  ’’Method  of 
random  testing  (the  Monte-Carlo  method)’’  transaction  of  the  Physics 
Inst,  of  the  Acad,  of  Sc.  Georgian  SSR,  Vol.  3,  pp.  107-122, 
"Statistical-probable  modeling  of  physical  processes  and  structures,’’ 
transactions  of  the  Inst,  of  Cybernetics  of  the  Acad,  of  Sc.  Georgian 
SSR,  Vol.  1,  pp.  13-18. 
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function  x(t),  containing  no  phase  information  and  bound  with  the 
spectral  density  of  the  prower  of  the  process  S(f)  by  the  Fourier 

transformation  '!  (*)  S  (j)e1  dj. 

-  -  a 

At  a  statistical-probability  modeling  of  any  random  process  the 
first  its  section,  equal  to  the  correlation  interval,1  is  represented 
in  form  of  a  system  of  n  random  values,  distributed  in  accordance  to 
the  given  monodimensional  function  F-,(x)  and  statistically  bound 
between  themselves  in  conformity  with  its  correlation  function  <J/(t)  at 
T0 

n  =  -££,  where  At  —  step  between  discrete  values  of  the  argument, 
corresponding  to  modeled  sections  of  the  process.  All  subsequent 
instantaneous  realization  values  of  the  process  from  moment 
t  =  t0  +  At  and  further  are  found  with  the  aid  of  a  drawing  of  random 
values  by  the  rule,  which  appears  to  be  an  n-dimensional  conditional 
distribution  function  and  is  determined  by  given  F^x)  and  *JKt)  at 
known  n  of  previous  instantaneous  values  of  that  realization. 

Considering  the  discrete  nature  of  statistical-probability 
modeling  of  random  processes  and  the  greater  volume  of  calculations, 
bound  with  the  formation  of  a  sequence  of  random  values  —  instantaneous 
realization  values,  for  mathematical  modeling  is  usually  used  an 
universal  or  special  electron  digital  computer.2 

The  below  discussed  two  instances  illustrated  practical  methods 
of  applying  the  statistical-probability  modeling  method. 


^he  correlation  interval  of  the  process  is  called  the  section  of 
time  Tq,  beyond  the  limits  of  which  statistical  bonds  between  sections 

of  the  process  are  practically  inexistent. 

zWhen  studying  noncomplex  situations,  when  the  drawing  is  made 
with  the  aid  of  random  number  tables,  the  numerical  calculation  can 
be  carried  out  also  with  the  aid  of  a  key  computer. 
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Calculation  of  Parameters  of  the  Process  of  Approaching 

of  a  Rocket  and  Target 

It  is  known,1  that  if  data,  obtained  from  the  system  of  target 
indication  to  assure  rocket  —  target  encounter,  are  insufficiently 
accurate,  and  the  rocket  has  a  highly  exact,  but  possessing  a  small 
radius  of  action  control  apparatus,  then  from  the  viewpoint  of 
calculated  evaluation  of  tactical  effectiveness  of  the  armament  a 
decisive  value  acquire  such  kinematic  parameters,  as  the  distance 
between  the  rocket  and  target,  bearing  to  the  target,  etc. 

We  will  calculate  the  minimum  distance  between  rocket  and  target 
for  the  following  conditions: 

the  target  moves  rectilinear ly  at  constant  velocity  Vfc; 

velocity  of  the  rocket  Vpi  for  every  i-realization  is  unchanged; 

movement  of  rocket  and  target  occurs  in  a  stationary  plane, 
formed  by  vectors  and 


the  moment  of  firing  (at  t 


The  relative  position  of  rocket  and 
target  at  each  moment  of  time  (t)  will  be 
determined  by  the  distance  between  same 
r(t),  angles  between  vectors  and  ^p 
and  by  the  direction  of  rocket  —  target. 

The  value  of  the  leading  angle  at 
0)  is  found  from  the  correlation 


IV 


(1) 


where  V. 

i 

target;  Vp 


and  q  —  measured  velocity  values  and  course  angle  of  the 
ui 

-  nominal  value  of  rocket  velocity. 


‘See  for  example,  A.  S.  Lokk.  Missile  control,  GITTL,  1957, 
pp.  327  and  498. 
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Results  of  measuring  target  movement  elements,  as  already 
mentioned,  appear  to  be  random  values,  which  are  distributed  by  the 
normal  law 


V„  *  K,,,  ou)  and  (}Uj  *  N (<7 «• 


where  Vfc  and  qfc  —  actual  values  of  important  magnitudes;  o^  and 
—  mean  quadratic  velocity  deviations  and  of  course  angle  from  average 


V. 


‘t* 


At  a  normal  distribution  of  random  value  X  with  a  probability, 
close  to  unity,  is  valid  ratio 


*5 


(2) 


where  AX _ _  -  maximum  possible  deviation  of  random  value  from  mean 

max 

value  x;  ax  -  corresponding  standard  or  mean  square  deviation. 


The  actual  value  of  rocket  velocity  V  and  angle  between  vector 

pi 

V  and  target  direction  at  an  i-discharge  differ  from  calculated  V  , 

<t>0,  because  they  appear  to  be  random  values,  which  are  distributed 

by  the  normal  law  V  N  (V  ,  a  )  and  <f>.  N  (<J>n,  a,),  where  a  and 

p*  p  i  0*  <p  p 

0^  —  mean  square  velocity  deviations  of  course  angle  of  the  rocket 

from  mean  V  and  4>n. 

P  U 


The  problem  about,  what  namely  random  value  magnitude  X  is  being 
realized  at  an  alternate  "theoretical  experiment,"  is  decided  either 
by  drawing  with  the  use  of  random  numbers  table,  or  with  the  aid  of 
a  sensing  element  of  pseudo  random  numbers. 


The  drawing  with  the  use  of  random  numbers1  is  based  on  a 


‘Detailed  information  about  this  is  given  in  the  above  mentioned 
report  of  V.  V.  Chavchanidze  "Method  of  random  tests." 


5 


prooability  ratio,  according  to  which  any  monodimensional  random  value 
X,  transformed  by  a  natural  integral  distribution  law  F(x),  becomes 
uniformly  distributed  in  the  range  of  from  0  to  1. 

When  using  tables  of  random  numbers  to  assure  independence  of 
selective  x.^  values  the  selection  of  random  numbers  for  all  situations 
is  necessary  to  realize  by  one  very  same  previously  established  law, 
beginning  with  any  arbitrary  selected  number. 

To  calculate  the  distance  between  rocket  and  target  at  any  moment 
of  time  (t)  it  is  possible  to  utilize  expression 

rt{i)  --  ]  Tv  -a jf-lW-f?.  ( 3 } 

where  rQ  —  distance  between  objects  at  moment  t  =  0, 

a, -  l'w  C03  Tn  v«  COS q,,.  b, ~  VPl  sin  «f,  \\  si n  qk. 

The  minimum  distance  between  rocket  and  target  at  i-discharge  is 
determined  from  conditions  =  0  and  is  expressed  by  formula 


_  _  f«]M_ 
mi 


(4) 


Carrying  out  a  sufficiently  great  number  of  "theoretical 
discharges"  M  and  having  determined  for  each  of  them  the  value  r^  it 

is  possible  to  calculate  the  statistical  probability  of  trapping  the 
target  at  any  given  radius  of  reaction  RQ  of  an  autonomous  control 
system1 


(5) 


1  Under  the  reaction  radius  RQ  is  understood  the  distance  of 
noncontact  detection  of  the  target  by  the  guidance  system. 
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where  m  —  number  of  "experiments,"  at  which  is  fulfilled  the  inequality 


rr,  4L 


Example.  Initial  data  are  given:  *  8  m/s,  ■  40°, 

Vp  =  25  m/s,  rQ  =  12,000  m,  AVt  =  1.8  m/s,  Aqt  *  18°,  AVp  *  0.6  m/s, 

\*i  W 

To  find  values  r_  for  the  given  situation. 

-  m± 

Solution  1.  We  determine  the  "measured"  values  of  elements  of 

target  movement  (V.  and  e.  ):  using  a  table  of  random  numbers,1  we 

zi  zi 

find,  that  0^  =  0806  and  02  =  5603;  with  the  aid  of  integral  functions 

table  for  the  normal  law  N  (0.1)  we  determine,  that  7X  »  *  —  i.w 

and  Z0  =  F**1  [0.56033  =  0.15; we  calculate  V.  ♦  N  (8;  0.6)  and  q. 

Z1  z± 

N  (40;  6)  and  we  find,  that  V.  =  0.6)  (-1.40)  +  8.0  =  7.2  m/s, 

ri 

q.  =  6.0-0.15  +  40  =  40.9°. 
zi 

2.  We  calculate  by  formula  (1)  value  of  the  lead  angle 


7,2  l 

•i0“arcsin - sin  {4C>,9  )[  ->0  50  . 

.  2a, 0  J 


3. 

rocket  V 

P 


We  determine  real  values  of  velocity  and  course  angle  of  the 

and  <j>.  for  which: 
i  1 


from  the  table  of  random  numbers  we  select  two  more  numbers,  we 
will  assume,  0^  a  9138  and  0^  =  0285;  by  table  of  integral  function 

we  determine  the  corresponding  values  Z?  =  F-^  (0.9138)  *  1.36  and 

—  1  * 

Zjj  =  F  ,  (0.0285)  =  -1.93;  we  calculate  values  of  sought  for 

parameters  of  rocket  motion 


-  t 

‘See  for  example,  A.  K.  Mitropol’ skiy .  —  Technique  of  statistical 
calculations.  Moscow,  Fizmatgiz,  1961. 
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I'p.  0.2  1.36  +  25,0  =  25,272  ■/*, 

H.  «  3  (  1.93)  +  IO'o>/— 5,0°. 

4.  We  calculate  by  formula  (4)  the  value  of  minimum  distance 
from  the  rocket  to  the  target  at  an  i-discharge 


12  003-2, frlf» 
«w67?5 


Reproducing  the  Realization  of  Noise  Process  with  Normal 

Law  of  Distribution 

The  physical  fields  spreading  in  real  media,  as  a  rule,  are  formed 
under  the  effect  of  many  factors,  having  a  random  nature.  As  result 
of  this  also  the  fields  represent  random  processes,  the  distribution 
of  which  is  subject  to  the  Gauss  law. 

At  statistical-probability  modeling  of  a  section  of  a  stationary 
random  process  with  normal  distribution  law  is  possible  to  use  the 
method  of  canonical  disposition  of  functions.1  Considering,  that 
really  existing  noise  processes  have  a  zero  mean  value,  canonical 
decomposition  of  corresponding  random  function  can  be  written  in  form 


x (/)-•■  V B,- \\0).  (6) 

K 

where  B,  —  random  noncorrelated  values  (decomposition  coefficients), 
distributed  by  normal  law  N  (0,  a^);  v^Ct)  “  coordinate  time  functions, 
bound  with  by  ratio 


V|(0“ 


(')J- 


(7) 


‘Detailed  data  about  this,  for  example,  in  the  V.  S.  Pugachev  book 
"Theory  of  Random  Functions,"  GIFML,  1962. 


FTD-HT-2  3- 101-68 


8 


If  is  presented  in  form  of  a  linear  combination  of  Instantaneous 
values  of  functions  x(t)  In  the  corresponding  points  of  the  modeled 
section  t0  and  formula  (7)  Is  considered,  then  is  possible  to  obtain,, 
recurrent  ratios  for  the  calculation  of  dispersions  of  coefficients 
Bt  and  coordinate  functions  vA(t) 

«t  l  ■  1  oi  v,(/)=r 

0, 

i  i 

at  t,)  -  1  oiKf/,))*  . 

«- 1 

vi(0"  A-  |  *,«.  /,)  -  *V  ,,i  .  vt(/)v,(/()j,  (8) 

.  ■  i 

where  ^  (t^,  t2>  —  correlation  function  of  the  process. 

i 

Practical  application  of  the  described  method  of  modeling  ‘  \ 
instantaneous  values  of  realising  a  random  function  will  be  examined} 

t  i 

0  » 
for  ^ -  =  3,  where  At  -  value  of  the  step  on  the  axis  of  time  between  . 

discrete  values  of  the  argument  in  section  points  of  the  process.  Af/ 

%  • 

first  are  made  auxiliary  calculation:;:  *  ‘l 

» 

1.  Calculated  are  values  of  elements  of  a  triangular  correlation 

matrix  |  tj>  (t)I  =  *».. v  <  ■  u)  | 

x  ‘  !,(")•  *,<•*'»  i 

i 

1 

2.  By  recurrent  ratios  (8)  are  calculated  three  values  o.  and 

I  * 

the  corresponding  to  the  correlation  matrix  number  of  elements  of 
triangular  matrix  ||v£(t)||(i  =  1,  2,  3). 

The  obtained  values  are  used  during  the  calculation  of  three 
instantaneous  values  of  function  x(t). 

1.  By  the  selected  method  (with  the  aid  of  table  of  random 
numbers  or  a  special  generator)  are  obtained  three  groups  of  random 
numbers  (12  numbers  in  every  group),  uniformly  distributed  in  the 
interval  (0,  1). 
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2.  Calculated  are  three  random  numbers  y,  distributed  by  normal 


law  JV(0,1)— Yi-Sp;  -6  (at  /-l, 2.3). 

3.  Determined  are  values  of  canonical  decomposition  coefficients 

B1  *  al*Yl»  B2  *  a2*Y2»  B3  “  a3*Y3*  where  ai»  a2»  a3  “  previously 
calculated  standard  deviations  of  random  values  of  coefficients 

B£  at  ' 

4.  During  the  multiplication  of  elements  of  a  triangular  matrix 
ii vi  (C  *i)!:  into  corresponding  values  is  obtained  a  triangular  matrix 

i!v,(/k.  <t) ■**!!. 

5.  Instantaneous  values  of  realizing  function  x(t)  on  section 
tq  is  obtained  by  adding  by  columns  of  elements  of  the  triangular 
matrix 


The  described  mathematical  method  allows  with  the  aid  of  a  TSVM1 
to  obtain  statistical  evaluations  of  characteristics  of  needed  accuracy 
and  to  reproduce  probability  models  of  investigated  objects  on  the 
basis  of  experimental  data  about  laws  of  distribution  and  correlation 
bonds  in  the  elements  of  these  objects,  not  expecting  the  formulation 
of  an  analytical  theory  and  revealing  their  physical  nature. 


^lnce  the  number  of  investigated  realizations  is  limited  only  by 
rapid  action  and  ’’memory”  volume  of  the  machine. 


’v* 
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